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Abstract 

We study finite nuclei, at the mean-field level, using the Zimanyi-Moskowski 
model and one of its variations (the ZM3 model). We calculate energy 
levels and ground-state properties in nuclei where the mean-field approach 
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is reliable. The role played by the spin-orbit potential in sorting out mean- 
field model descriptions is emphasized. 
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I. INTRODUCTION 



At the mean-field level the linear o — uj (or Walecka) model Q satisfactorily explains 
many properties of nuclear matter and finite nuclei with two free parameters. The re- 
sulting nuclear-matter compression modulus at saturation density exceeds, however, the 
experimental bound A way out of this difficulty is to introduce non-linear scalar 

self-couplings |||| . The non-linear model that obtains has been shown to reproduce well 
ground-state nuclear properties (though with some instabilities at high densities and for 
low values of the compressibility, k < 200 MeV BJ). This model is renormalizable and 
has four free parameters to fit. 

One alternative approach which renders a satisfactory compression modulus without 
increasing the number of free parameters is that advanced by Zimanyi and Moszkowski 
and by Heide and Rudaz 0. This model employs the same degrees of freedom and the 
same number of independent couplings that are present in the Walecka model. The 
difference lies in the coupling among the fields. In the work of Zimanyi and Moszkowski 
(ZM), for instance, the authors use a non-renormalizable derivative coupling between the 
scalar-meson and the baryon fields which they later adjust to reproduce the experimental 
conditions at saturation. Their results for the compression modulus, k = 224 MeV, and 
effective mass, M* = 797 MeV, compare very well with Skyrme-type calculations ||. 

Since it is desirable to have models that, at the mean-field level, provide reasonable 
nuclear matter results, the original ZM model and its two variations (to which we shall 
refer hereafter as ZM2 and ZM3) described in the appendix of the original paper of 
Zimanyi and Moszkowski M, have deserved some recent attention. Particularly, Koepf, 
Sharma and Ring || (KSR) have performed nuclear-matter and finite-nuclei calculations 
and compared the ZM model against i) the linear a — uj model and, ii) models with non- 
linear scalar couplings of phenomenological origin. One of their findings was that the ZM 
model, when used to calculate the energy spectrum of 208 Pb, gives a spin-orbit splitting 
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that compares poorly with the data. It is well known that the splitting of the energy levels 
due to the spin-orbit interaction is very sensitive to the difference X=V — S between the 
vector, V, and the scalar, S, potentials. For the ZM model one gets A=223 MeV, which 
is short of the 785 MeV of the successful (in this regard) Walecka model. A simplified 
summary of KSR's results also indicates that: 1) The (non)linear a — uo model predicts 
a (good)reasonable spin-orbit splitting, a (soft)stiff equation of state and a (good)small 
effective nucleon mass; 2) The ZM model predicts a soft equation of state, a reasonable 
effective nucleon mass and a poor spin-orbit splitting for finite nuclei. To this list we may 
add [ |TU| , |rTf that the relativistic content of the various models, as given by the ratio of 
scalar to baryon densities, differs also being the linear and non-linear o — uo models more 
relativistic than the usual ZM. 

At this point it is clear that fixing the spin-orbit problem in the ZM model would 
make it similar in results to the non-linear a — u model but with the additional bonus of 
having to fit two parameters instead of four. The modified ZM models, ZM2 and ZM3, 
aim at this. All three models come about in the following way. In the standard ZM model 
the non-linear effective scalar factor l/m*=l + g a a/M, whose meaning will become clear 
below, multiplies the nucleon derivative and the nucleon- vector coupling terms. When its 
two free parameters are fitted to the nuclear-matter baryon density, po=0.148 fm -3 , and 
energy density per nucleon at saturation, E^= — 15.75 MeV, one obtains the results shown 
in the first row of Table |. If we now let 1 / m* to act upon the nucleon derivative and all 
terms involving the vector field we end up with the ZM2 version of the ZM model, whose 
predictions are shown in the second row of the table. The compressibility is slightly lower 
and the difference between the vector and scalar potentials has increased to 278 MeV [l I 



Finally, the ZM3 version of the model is obtained by allowing 1/m* to act just on the 
nucleon derivative term. With this change one obtains the results displayed in the third 



row of the table; the compressibility diminishes again and A reaches 471 MeV [11]. It is 



4 



worth calling the attention to the fact that the particular non-linear couplings used in these 
models change the traditional conception whereby to a lower value of m* = M* /M should 
correspond a higher compressibility k. Particularly, the ZM3 field equations couple the 
non-linear scalar field to the vector field and this new source is responsible for the m*=0.72 
and k=156 MeV, of the table. In the phenomenological non-linear scalar-coupling model 
of Feldmeir and Lindner ]I2[ (hiperbolic ansatz), for example, the price to pay for the 
reduced effective mass, m* = 0.71, is an enhanced compressibility, k = 410 MeV. 

If we assume that the difference A between the scalar and the vector potentials is a 
good indicator of the strength of the spin-orbit splitting in finite nuclei, it is clear from 
Table | that ZM3 offers the closest approach of the three to the non-linear a — uo model. 
In this situation it is relevant to study the model predictions for a number of properties 
of finite nuclei for which there exist experimental data. In this work we apply the ZM3 
model to the study of the energy spectrum and ground-state properties of 16 0, 40 Ca, 
48 Ca, 90 Zr and 208 Pb. We calculate energy levels, charge r.m.s radius and nucleon density 
distributions and compare the results obtained with the predictions of the standard ZM 
model and with the experimental data in order to extract conclusions. 

In the next section we introduce the ZM and ZM3 models and give the necessary detail 
to understand the origin of the results we obtain. These are presented in the last section 
where we also draw some conclusions. 



II. THE MODEL 



The ZM and ZM3 models of interest for this work can be derived from the Lagrangean 



a=l 



id" 



(M-^»Ua (1) 
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where the effective coupling constants are given in each model according to 



Model g* 




ZM m*g a 




ZM3 m*g a 


m*g^p 



where we define 



/ i i ' ' 

m = 1 + 



M 

In equation ([[J) the nucleon field, ip, couples to the scalar-isoscalar meson field, a, and 
the vector-isoscalar meson field, uj^. A third isovector meson field, p M , (neutral component) 
is included to account for the asymmetry between protons and neutrons. The p M and the 
electromagnetic field both couple to the baryon field. 

The Euler-Lagrange equations of motion using the Lagrangean (|I|) give the following 
equations for the fields, 



In 



-(M-0»Ua = O, (2) 



dyvT + rrfo v = gl E ^7>a , (3) 

a=l 

d»(T + m 2 ^ = ^Y1 ^al V T^ a , (4) 
Z 0=1 

= |E + , (5) 

^ a=l 

A A /do* - Ida* - \ 

d^a + m\u = m* 2 g a E Ma - E ( ffiaJ^aU 1 * + 2 ■ffiaJvTtfpaff) • (6) 

In the mean-field approximation all baryon currents are replaced by their ground- 
state expectation values. In a system with spherical symmetry the mean value of the 
spatial components of the vector-meson fields vanish, resulting in the following mean-field 
equations, 
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Equations (0) to (p"T|) are a set of coupled non-lineal differential equations which may 
be solved by iteration. For a given set of initial meson potentials, the Dirac equation (|7|) 
is solved. Once the baryon wave functions are determined the source terms for Eqs. (|8|)- 
(|TT|) are evaluated and new meson fields are obtained. The procedure is iterated until 
self-consistency is achieved. 

The coupling constants g„ and g^ are chosen to reproduce the saturation baryon 
density in symmetric nuclear matter, po=0.148 fm~ 3 , and the energy density per nucleon 
at saturation, E/A— — 15.75 MeV. The third coupling constant, g p , is obtained by fitting 
the bulk symmetry energy in nuclear matter given by the expression 



9 2 P (g;\ kl 



im P \9p 

to a 4 =32.5 MeV at p=po- We have used 



(12) 



Po 



Q 2 F 



and 



E* = Jkl + m* 2 M 2 . 



The coupling constants and masses used in the calculations are given in Table |I|. 



III. RESULTS 



Results of the present calculation for the single particle energy levels in 16 0, 40 Ca, 
48 Ca, 90 Zr, and 208 Pb are shown in Tables JT| to [V] Through the tables a good overall 
agreement is obtained, i) between the two models and, ii) with the experimental data. 
Nonetheless, the ZM3 model fares better than ZM in reproducing the observed spin-orbit 
splittings. This is due to the difference in the way that the scalar and vector mesonic 
fields couple to each other — via non-linear effective coupling constants — in each model. 
This difference is better illustrated in figures |l] to § where the "central potential", Vo, and 
the "reduced spin-orbit term", Vi s , are depicted as functions of the radial distance, r, for 
oxygen and lead. We have defined the central potential to be 

V = V + S (13) 

with 

V = g*wo 

and 

S = -g>- 

The reduced spin-orbit term, on the other hand, is given by the expression, 

V - S' 

Vis = - , (14) 

2m s r 

with m s = M — 0.5* (V — S) and with the primes denoting partial derivatives with respect 
to r. Though the choice of mass in the denominator of Eq. ([14]) could be opened to some 
debate, the figures are intended to compare strengths in both models and not to draw 
absolute conclusions. 

Figures [l] and ^| illustrate the origin of the behaviour of both models. ZM is shallower 
in the center and raises more steeply on the surface than ZM3; this holds true both in 
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oxygen and in lead. The resulting spin-orbit potential (figures |] and £|) is, thus, deeper at 

the surface in ZM3 than in ZM. To estimate the magnitude of this difference we calculated 

the ratio expected for the spin-orbit splittings in ZM3 and ZM, using the parameters from 

the fit to nuclear matter of Table |. The result 

Vi Z s M " (V-S) z ^m* ZM 
v zm <* {v _ S )^ m * ZM3 ~ ^ 

agrees nicely with the ratios extracted from the levels in tables |ITT|-|V|. 

It is worth noticing that our calculated spectra for the ZM model, differ from those 
presented by KSR in Ref. || for the same model. We have traced the origin of this 
discrepancy to the symmetry energy used to fit the p-meson coupling constant. Figure [5] 
illustrates the effect that a decreasing 04 has on the energy of a few selected single particle 
orbits in 208 Pb. The shift in energy due to the presence of the p-meson may be as large as 5 
MeV depending on the choice of 04 (the arrow indicates the value used in the calculations 
shown here, 04 = 32.5 MeV). The spectrum of 208 Pb calculated with a vanishing g p agrees 
with the results of KSR. Notice also that the splittings amongst spin-orbit partners are 
insensitive to g p since changes in 04 shift the levels globally. 

At this point one of our conclusions is that, as demonstrated in ||, the ZM model 
does not do well in reproducing the energy splittings due to the spin-orbit interaction. 
However, the overall spectrum turns out to be satisfactory. In this sense the isovector 
meson plays in ZM the same important role, for asymmetric nuclei, that has been shown 
to play in the non- linear Walecka model fi"3||. 



Table |\T| shows the results obtained for some of the static ground-state properties of 
the same nuclei as above. One feature to point out is that, systematically, ZM3 predicts 
a r.m.s. value for the charge radius that is larger than the one calculated with the ZM 
model. This is due to the slightly smaller binding of the protons in this model that 
produces a longer tail in the charge distribution, as shown in Figs. § and [7|. The baryon 
density extends also farther in ZM3 as can be appreciated from Figs. [8] and || The baryon 



distribution at the surface indicates that edge effects are more important in ZM3 than in 
the ZM model, something that agrees with our discussion of the spin-orbit splittings in 
the previous paragraphs. 

Summarizing, we have applied the ZM3 model to the study of the energy spectrum and 
ground-state properties of 16 0, 40 Ca, 48 Ca, 90 Zr and 208 Pb. The interest in this model, and 
the standard ZM, resides in their ability to describe nuclear matter at saturation with two 
free parameters. For finite nuclei, we have shown that ZM3 gives a reasonable description 
of nuclear spectra and improves upon ZM regarding the energy splitting of levels due to 
the spin-orbit interaction. The results of the calculations described here were compared 
with those obtained using the standard ZM model and with the experimental data. It is 
our conclusion that models of this kind, with derivative couplings involving the scalar and 
the vector fields, offer a valid framework to pursue calculations where the requirement of 
simultaneous reasonable values for the nuclear compressibility and the spin-orbit splitting 
cannot be side-stepped. 
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TABLES 

TABLE I. Nuclear matter results in the three Zimanyi-Moskowszki models. The compress- 
ibility k, the scalar potential S, and the vector potential V, are in MeV. The effective mass is 

in units of the bare mass. 

Model k S V m* 

ZM 225. -141. 82.5 0.85 

ZM2 198. -168. 110.7 0.82 

ZM3 156. -267. 204.7 0.72 



TABLE II. Parameters used in the calculation. 



Model 


m N 


m a 


m w 


m p 


9l 


al 


9l 


ZM 


938.27 


520 


783 


770 


55.540 


44.207 


79.605 


ZM3 


938.27 


520 


783 


770 


135.30 


211.21 


137.74 



12 



TABLE III. Single-particle energy levels (in MeV) calculated in the Zimanyi-Moszkowski 
model (ZM) and in the ZM3 version. Where available the experimental energies are also quoted. 
Numbers between brackets correspond to protons, others to neutrons. 



Nucl. Levels ZM ZM3 Exp. 



16Q 


ls l/2 


35.4 


(31.2) 


36.2 


(32.1) 


47.0 


(40±8) 




1P3/2 


19.6 


(15.6) 


19.4 


(15.6) 


21.8 


(18.4) 




lPl/2 


18.2 


(14.2) 


16.5 


(12.7) 


15.7 


(12.1) 


40 Ca 


lSl/2 


43.3 


(35.1) 


46.7 


(38.7) 




(50±10) 




103/2 


32.4 


(24.5) 


33.6 


(25.9) 




(34±6) 




lPl/2 


31.6 


(23.7) 


31.9 


(24.2) 




(34±6) 




1^5/2 


20.2 


(12.6) 


20.4 


(13.0) 


21.9 


(15.5) 




1^3/2 


18.7 


(11.2) 


17.3 


(9.97) 


15.6 


(8.3) 




2si/ 2 


15.9 


(8.33) 


16.9 


(9.52) 


18.2 


(10.9) 


48 Ca 


lSl/2 


41.8 


(39.9) 


45.7 


(43.2) 




(55±9) 




1P3/2 


31.3 


(30.7) 


33.2 


(32.1) 




(35±7) 




lPl/2 


30.8 


(30.1) 


32.0 


(30.7) 




(35±7) 




1^5/2 


19.8 


(19.7) 


20.4 


(20.0) 


16.0 


(20.0) 




1^3/2 


18.7 


(18.5) 


17.9 


(17.4) 


12.4 


(15.3) 




1/7/2 


7.87 




7.94 




9.9 






2si/ 2 


16.7 


(14.4) 


17.2 


(15.4) 


12.4 


(15.8) 
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TABLE IV. Same as Table 



HI. 



Nucl. 


Levels 




ZM 




ZM3 




Exp. 


90 Zr 


Lsi/2 


45.9 


(37.2) 


51.0 


(41.6) 




(54±8) 




J-P3/2 


38.4 


(30.8) 


41.6 


(33.5) 




(43±8) 




lj»l/2 


38.1 


(30.5) 


40.9 


(32.7) 




(43±8) 




14/2 


29.7 


(22.9) 


31.4 


(24.3) 




(27±8) 




14/2 


29.1 


(22.2) 


29.9 


(22.6) 




(27±8) 




1/7/2 


20.3 


(13.8) 


20.8 


(14.3) 








1/5/2 


19.1 


(12.6) 


18.3 


(11.7) 


13.5 






109/2 


10.2 




10.2 




12.0 






2si/ 2 


26.6 


(18.2) 


28.1 


(19.9) 




(27±8) 




2p 3 /2 


16.0 


(7.64) 


16.9 


(8.88) 


13.1 






2P1/2 


15.6 


(7.25) 


16.0 


(8.01) 


12.6 
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TABLE V. Same as Table 



ID. 



Nucl. 


Levels 




ZM 




ZM3 


Exp. 


208p b 




45.1 


(35.6) 


50.6 


(40.0) 






J-P3/2 


41.0 


(31.8) 


45.3 


(35.2) 






lPl/2 


40.7 


(31.6) 


45.1 


(34.9) 






14/2 


35.8 


(26.9) 


39.0 


(29.4) 






14/2 


35.6 


(26.6) 


38.4 


(28.7) 






1/7/2 


29.9 


(21.2) 


31.9 


(22.7) 






1/5/2 


29.4 


(20.6) 


30.9 


(21.5) 






159/2 


23.2 


(14.6) 


24.3 


(15.4) 






157/2 


22.4 


(13.8) 


22.6 


(13.6) 


(11.4) 




1^11/2 


15.9 


(7.41) 


16.2 


(7.65) 


(9.4) 




l/lg/2 


14.8 




13.9 




10.8 




1A3/2 


8.18 




7.98 




9.0 




2^1/2 


33.2 


(23.7) 


36.2 


(26.1) 






2p 3 /2 


25.8 


(16.5) 


27.8 


(18.1) 






2pi/ 2 


25.6 


(16.3) 


27.4 


(17.6) 






24/2 


18.1 


(8.87) 


19.4 


(9.92) 


(9.7) 




24/2 


17.7 


(8.50) 


18.5 


(9.06) 


(8.4) 




2/7/2 


10.2 




11.0 




9.7 




2/5/2 


9.69 




9.85 




7.9 




3si/ 2 


16.3 


(6.75) 


17.6 


(7.81) 


(8.0) 




3p 3 /2 


7.88 




8.83 




8.3 




3pi/ 2 


7.69 




8.38 




7.4 
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TABLE VI. The binding energy per nucleon e (in MeV), the mean square charge radius (r 2 ) 
(in fm 2 ) and, the spin-orbit splitting (for partners where data are available) AE so (MeV), are 
shown for the ZM and ZM3 models. The experimental values are shown for reference. 



Nucleus g.s.p. ZM ZM3 Exp. 

16 e -8.40 -7.50 -7.98 

(r 2 ) 2.64 2.78 2.74 

AE so (lp 3/2 -lp 1/2 ) 1.4 2.8 6.1 

40 Ca e -9.01 -8.43 -8.55 

(r 2 ) 3.39 3.51 3.48 

AE so (ld 5/2 -ld 3/2 ) 1.5 6.1 6.3 

48 Ca £ -8.78 -8.30 -8.67 

(r 2 ) 3.47 3.57 3.47 

AE so (ld 5/2 -ld 3/2 ) 1.1 2.5 3.6 

90 Zr £ -8.80 -8.46 -8.71 

(r 2 ) 4.25 4.35 4.27 

AE so (2p 3/2 -2 Pl/2 ) 0.4 0.9 0.5 

208 Pb £ -7.86 -7.66 -7.87 

(r 2 ) 5.54 5.66 5.50 

AE so (2f 7/2 -2f 5/2 ) 0.5 1.15 1.8 

AE so (3p 3/2 -3 Pl/2 ) 0.19 0.45 0.9 
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FIGURES 

FIG. 1. The potential Vq (denned in equation ( |l3|) in the text) as a function of the ra- 
dial distance. The calculation was done in oxygen and the results correspond to the standard 
Zimanyi-Moszkowski model (ZM) Q and to the ZM3 model studied in this paper. 

FIG. 2. Same as figure p] but for the 208 Pb nucleus. 

FIG. 3. The spin-orbit potential V[ s (defined in equation (|l4|) in the text) as a function of the 
radial distance. The calculation was done in oxygen and the results correspond to the standard 
Zimanyi-Moszkowski model (ZM) Q and to the ZM3 model studied in this paper. 

FIG. 4. Same as figure || but for the 208 Pb nucleus. 

FIG. 5. Valence proton and neutron orbitals in 208 Pb as a function of the assymetry param- 
eter 04. The arrow indicates the value used in the calculations performed in this work. 

FIG. 6. The charge density distribution in 16 as a function of the radial distance. The 
results correspond to the standard Zimanyi-Moszkowski model (ZM) || and to the ZM3 model 
studied in this paper. 

FIG. 7. Same as figure ^ but for the 208 Pb nucleus. 

FIG. 8. The baryon density distribution in 16 as a function of the radial distance. The 
results correspond to the standard Zimanyi-Moszkowski model (ZM) || and to the ZM3 model 
studied in this paper. 



FIG. 9. Same as figure || but for the 208 Pb nucleus. 
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